Fractional ac Josephson effect in unconventional superconductors 
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For certain orientations of Josephson junctions between two p x -wave or two d-wave supercon- 
ductors, the subgap Andreev bound states produce a 47r-periodic relation between the Josephson 
current I and the phase difference <f>: I oc sin(0/2). Consequently, the ac Josephson current has the 
fractional frequency eV/h, where V is the dc voltage. In the tunneling limit, the Josephson current is 
proportional to the first power (not square) of the electron tunneling amplitude. Thus, the Joseph- 
son current between unconventional superconductors is carried by single electrons, rather than by 
Cooper pairs. The fractional ac Josephson effect can be observed experimentally by measuring 
frequency spectrum of microwave radiation from the junction. 

PACS numbers: 74.50.+r 74.70.Kn 74.72.-h 



I. BRIEF HISTORY OF THE AC JOSEPHSON 
EFFECT 

In 1962, Josephson [1] predicted theoretically that if a 
dc voltage V is applied to a junction between two super- 
conductors, ac supercurrent with the frequency 2eV/h 
appears between the superconductors. The ac Josephson 
radiation was first observed experimentally 40 years ago 
in Kharkov by Yanson, Svistunov, and Dmitrenko [2, 3]. 
In Ref. [3], the spectrum of microwave radiation from 
tin junctions was measured, and a sharp peak in the fre- 
quency spectrum at 2eV/h was found. It is amazing that 
without any attempt to match impedances of the junc- 
tion and waveguide, Dmitrenko and Yanson [3] found the 
signal several hundred times stronger than the noise and 
the ratio of linewidth to the Josephson frequency less 
than 10~ 3 . This discovery was followed by further work 
in the United States [4] and Ukraine [5]. The results of 
these investigations were summarized in the classic book 
[6]. Since then, the ac Josephson radiation has been ob- 
served in many materials in various experimental setups. 
For example, a peak of Josephson radiation was found 
in Ref. [7] in indium junctions at the frequency 9 GHz 
with the width 36 MHz. In Ref. [8], a peak of Josephson 
radiation was observed around 11 GHz with the width 50 
MHz in Bi2Sr2CaCu208 single crystals with the current 
along the c axis perpendicular to the layers. 

The theory of the Josephson effect was originally de- 
veloped for conventional s-wavc superconductors. In this 
paper, we study Josephson junctions between unconven- 
tional superconductors, such as d-wave cuprates or p x - 
wave organic superconductors. We show that the midgap 
Andreev states in these materials produce a 47r-periodic 
relation between the Josephson current I and the phase 
difference 4>: I oc sin(</>/2). Consequently, the ac Joseph- 
son current has the fractional frequency eV/K, a half of 
the conventional value. We hope that this effect can ob- 
served experimentally as the corresponding peak in the 
frequency spectrum of Josephson radiation from uncon- 
ventional superconductors, such as d-wave cuprates, in 
the manner similar to the pioneering experiments [2, 3] 



performed on conventional s-wave superconductors. 



II. INTRODUCTION 

In many materials, the symmetry of the supercon- 
ducting order parameter is unconventional, i.e. not s- 
wave. In the high-T c cuprates, it is the singlet d x 2_ y 2- 
wave [9]. There is experimental evidence that, in the 
quasi-one-dimensional (Q1D) organic superconductors 
(TMTSF) 2 X [10], the symmetry is triplet [11], most likely 
the p^-wave [12], where the x axis is along the con- 
ducting chains. The unconventional pairing symmetry 
typically results in formation of midgap Andreev bound 
states on the surfaces of these superconductors. For 
d-wave cuprate superconductors, the midgap Andreev 
states were predicted theoretically in Ref. [13] and dis- 
covered experimentally as a zero-bias conductance peak 
in tunneling between normal metals and superconductors 
(see review [14]). For the Q1D organic superconductors, 
the midgap states were theoretically predicted to exist 
at the edges perpendicular to the chains [15, 16]. When 
two unconventional superconductors are joined together 
in a Josephson junction, their Andreev surface states hy- 
bridize to form Andreev bound states in the junction. 
These states are important for the Josephson current. 
Andreev bound states in high-T c junctions were reviewed 
in Ref. [17]. The Josephson effect between two Q1D p- 
wave superconductors was studied in Refs. [18, 19]. 

In the present paper, we predict a new effect for 
Josephson junctions between unconventional (nonchiral) 
superconductors, which we call the fractional ac Joseph- 
son effect. Suppose both superconductors forming a 
Josephson junction have surface midgap states originally. 
This is the case for the butt-to-butt junction between two 
p^-wave Q1D superconductors, as shown in Fig. la, and 
for the 45°/45° in-plane junction between two d-wave 
superconductors, as shown in Fig. 3a. (The two angles 
indicate the orientation of the junction line relative to 
the b axes of each d x 2_ y a superconductor.) We predict 
that the contribution of the hybridized Andreev bound 
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states produces a 47r-periodic relation between the su- 
percurrent I and the superconducting phase difference 
(j>: I oc sm(4>/2) [20]. Consequently, the ac Josephson 
effect has the frequency eV/U, where e is the electron 
charge, V is the applied dc voltage, and K is the Planck 
constant. The predicted frequency is a half of the conven- 
tional Josephson frequency 2eV/h originating from the 
conventional Josephson relation I oc sin <j> with the period 
of 2ir. Qualitatively, the predicted effect can be under- 
stood as follows. The Josephson current across the two 
unconventional superconductors is carried by tunneling 
of single electrons (rather than Cooper pairs) between 
the two resonant midgap states. Thus, the Cooper pair 
charge 2e is replaced the single charge e in the expression 
for the Josephson frequency. This interpretation is also 
supported by the finding that, in the tunneling limit, the 
Josephson current is proportional to the first power (not 
square) of the electron tunneling amplitude [21-23]. Pos- 
sibilities for experimental observation of the fractional ac 
Josephson effect are discussed in Sec. VI. 

The predicted current-phase relation I oc sin(c6/2) is 
quite radical, because every textbook on superconductiv- 
ity says that the Josephson current must be a 27r-periodic 
function of <f> [20]. To our knowledge, the only paper that 
discussed the 47r-periodic Josephson effect is Ref. [24] by 
Kitaev. He considered a highly idealized model of spin- 
less fermions on a one-dimensional (ID) lattice with su- 
perconducting pairing on the neighboring sites. The pair- 
ing potential in this case has the p x -wave symmetry, and 
midgap states exist at the ends of the chain. They are 
described by the Majorana fermions, which Kitaev pro- 
posed to use for nonvolatile memory in quantum comput- 
ing. He found that, when two such superconductors are 
brought in contact, the system is 47r-periodic in the phase 
difference between the superconductors. Our results are 
in agreement with his work. However, we formulate the 
problem as an experimentally realistic Josephson effect 
between known superconducting materials. 



III. THE BASICS 

In this paper, we consider singlet pairing and triplet 
pairing with the spin polarization vector n having a uni- 
form, momentum-independent orientation [11, 12]. If the 
spin quantization axis z is selected along n, then the 
Cooper pairing takes place between electrons with the op- 
posite z-axis spin projections a and a: {c a {k)c s (—k)) oc 
A (J (fc), where c a (k) is the annihilation operator of an 
electron with momentum k and spin a. The pairing po- 
tential has the symmetry A CT (fc) = ^fAg-(k) — ±A a (— k), 
where the upper and lower signs correspond to the singlet 
and triplet cases. 

We select the coordinate axis x perpendicular to the 
Josephson junction plane. We assume that the interface 
between the two superconductors is smooth enough, so 
that the electron momentum component k y parallel to 
the junction plane is a conserved good quantum number. 



Electron states in a superconductor are described by 
the Bogoliubov operators 7, which are related to the elec- 
tron operators c by the following equations [25] 

c-k y (a) - (x) + v n*ky ^ (2) 

n 

where k y — —k y , and n is the quantum number of 
the Bogoliubov eigenstates. The two-components vec- 
tors ipnak y (x) = [u n(yky (x),v naky (x)\ are the eigenstates 
of the Bogoliubov-de Gennes (BdG) equation with the 
eigenenergies E naky : 

(e k (k x ) + U(x) A ak {x,k x ) \ ._ E . (3) 
{ Al ky (x,k x ) - £ky (k x )-U(x) [6 > 

where k x = —id x is the x component of the electron 
momentum operator, and U(x) is a potential. In Eq. (3) 
and below, we often omit the indices a and k y to shorten 
notation where it does not cause confusion. 



IV. JUNCTIONS BETWEEN 
QUASI-ONE-DIMENSIONAL 
SUPERCONDUCTORS 

In this section, we consider junctions between two 
Q1D superconductors, such as organic superconductors 
(TMTSF)2X, with the chains along the x axis, as shown 
in Fig. la. For a Q1D conductor, the electron energy 
dispersion in Eq. (3) can be written as e = Ti 2 k 2 /2m — 
2tbCOs(bk y ) — /i, where m is an effective mass, /i is the 
chemical potential, b and tt are the distance and the tun- 
neling amplitude between the chains. The superconduct- 
ing pairing potentials in the s- and p^-wave cases have 
the forms 

A ( h \ / s-wave, 
A aky {x,k x ) - \ Ap k x/kFj px . wave; W 

where hkF = \/2m(i is the Fermi momentum, and a is 
treated as + for \ and — for J.. The index /3 = R,L labels 
the right (x > 0) and left (a; < 0) sides of the junction, 
and Af) acquires a phase difference <p across the junction: 

A L = A , A R = Aoe*. (5) 

The potential U(x) = UqS(x) in Eq. (3) represents the 
junction barrier located at x — 0. Integrating Eq. (3) 
over x from -0 to +0, we find the boundary conditions 
at x = 0: 

= ipR, d x tp R - d x ip L = k F Zip(0), (6) 
Z = 2mU a /h 2 k F , D = A/{Z 2 + 4), (7) 

where D is the transmission coefficient of the barrier. 
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A. Andreev bound states 

A general solution of Eq. (3) is a superposition of the 
terms with the momenta close to afcp, where the index 
a = ± labels the right- and left-moving electrons: 



(a) 



A 



Aoe" 



A 



U/3<T+ 



ikpx 



-ikpx 



(8) 

Here (3 = =F for R and L. Eq. (8) describes a subgap 
bound state with an energy \E\ < A , which is local- 
ized at the junction and decays exponentially in x within 
the length 1/k. The coefficients (w/3 CT a, «/3cra) in Eq. (8) 
arc determined by substituting the right- and left-moving 
terms separately into Eq. (3) for x ^ 0, where U(x) = 0. 
In the limit kp 3> n 1 we find 



U(3o 



Vfjo 



E + iafihKVF ' 



TiVF 



(9) 



where vf = fikp/m is the Fermi velocity, and Ap aa is 
equal to aAp for s-wave and to aAp for p x -wave, with A^ 
given by Eq. (5). The fc y -dependent Fermi momentum 
UUf — hkp + 2ti,cos(bky)/vF in Eq. (8) eliminates the 
dispersion in k y from the BdG equation. 

Substituting Eq. (8) into the boundary conditions (6), 
we obtain the linear homogeneous equations for the co- 
efficients A/3 and Bp. The compatibility condition for 
these equations gives an equation for the energies of the 
Andreev bound states. There are two subgap states with 
the energies E a — aE ((f>) labeled by the index a = ±: 



E { Q s) (6) = -A \Jl - L>sin 2 0/2), s-s junction,(10) 
Eo\4>) = -A \ r Dcos(<f>/2) 7 p x -p x junction. (11) 



The energies (10) and (11) are plotted as functions 
of 6 in the left panels (b) and (c) of Fig. 1. Without 
barrier (D = 1), the spectra of the s-s and p x -p x junc- 
tions arc the same and consist of two crossing curves 
E = =fAo cos(</>/2), shown by the thin lines in the left 
panel of Fig. lb. A non-zero barrier (D < 1) changes 
the energies of the Andreev bound states in the s-s and 
Px'Px junctions in different ways. In the s-s case, the two 
energy levels repel near (p — 7r and form two separated 
27r-periodic branches shown by the thick lines in the left 
panel of Fig. lb [25, 26]. In contrast, in the p x -p x case, 
the two energy levels continue to cross at <p — an d 
they separate from the continuum of states above +Ao 
and below — Ao, as show in the left panel of Fig. lc. The 
absence of energy levels repulsion indicates that there is 
no matrix element between these levels at <p — m the 
Px-Px case, unlike in the s-s case. 

As shown in Sec. V, the 45°/45° junction between two 
d-wave superconductors is mathematically equivalent to 
the p x -p x junction. Eq. (11) was derived for the 45°/45° 
junction in Ref. [22, 23, 27]. 
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FIG. 1: (a) Josephson junction between two Q1D p^-wave 
superconductors, (b) The energies (left panel) and the cur- 
rents (right panel) of the subgap states in the s-s junction 
as functions of the phase difference cj> for D = 1 (thin lines) 
and D = 0.9 (thick lines), (c) The same as (b) for the p x -p x 
junction at D — 0.2. 



B. dc Josephson effect in thermodynamic 
equilibrium 

It is well known [25, 28] that the current carried by a 
quasiparticle state a is 



2e dE a 



(12) 



The two subgap states carry opposite currents, which 
are plotted vs. <p m the right panels (b) and (c) of Fig. 
1 for the s-s and p x -p x junctions. In thermodynamic 
equilibrium, the total current is determined by the Fermi 
occupation numbers f a of the states at a temperature T: 



2e dE a 

n ^-f 86 u 



a=± 



¥^ tan H^j- (13) 



For the s-s junction, substituting Eq. (10) into Eq. (13), 
we recover the Ambegaokar-Baratoff formula [29] in the 
tunneling limit D <C 1 

Is „ D Sln _ tanh l_ j = Sln ^ __ tanh ^- ) 

(14) 

and the Kulik-Omelyanchuk formula [30] in the transpar- 
ent limit D — > 1 



sm 



^^tanhf AoCOS ^V 
2 J h V 2T J 



(15) 



Taking into account that the total current is proportional 
to the number TV of conducting channels in the junction 
(e.g. the number of chains), we have replaced the trans- 
mission coefficient D in Eq. (14) by the junction resis- 
tance R = h/2Ne 2 D in the normal state. 
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T/T c 



of (f>. The assumption of thermodynamic equilibrium was 
explicitly made in Rcf. [22] and was implicitly invoked 
in Refs. [18, 19, 21] by using the Matsubara diagram 
technique. In Ref. [31], temperature dependence of the 
Josephson critical current was measured in the YBCO 
ramp-edge junctions with different crystal angles and was 
found to be qualitatively consistent with the upper curve 
in Fig. 2. 



C. Dynamical fractional ac Josephson effect 



FIG. 2: Critical currents of the s-s (dashed line) and p x -p x 
(solid line) Josephson junctions as functions of temperature 
for D = 0.3. 

Substituting Eq. (11) into Eq. (13), we find the Joseph- 
son current in the p x -p x junction in thermodynamic equi- 
librium: 

The temperature dependences of the critical currents 
for the s-s and p x -p x junctions are shown in Fig. 2. They 
are obtained from Eqs. (14) and (16) assuming the BCS 
temperature dependence for Ao. In the vicinity of T c , 
I p and I s have the same behavior. With the decrease 
of temperature, I s quickly saturates to a constant value, 
because, for D < 1, E ( a s) w =fA (10), thus, for T < A , 
the upper subgap state is empty and the lower one is com- 
pletely filled. In contrast, I p rapidly increases with de- 
creasing temperature as 1 /T and saturates to a value en- 
hanced by the factor 2/VZ) relative to the Ambegaokar- 
Baratoff formula (10) at T = 0. This is a consequence 
of two effects. As Eqs. (14) and (16) show, I s oc D and 
I p oc >/D, thus I p ^> I s in the tunneling limit D<1. At 
the same time, the energy splitting between the two sub- 
gap states in the p x -p x junction is small compared to the 
gap: £<j p) oc VdAq « A . Thus, for V~DA < T < A , 
the two subgap states are almost equally populated, so 
the critical current has the 1/T temperature dependence 
analogous to the Curie spin susceptibility. 

Eq. (16) was derived analytically for the 45°/45° 
junction between two c?-wave superconductors in Refs. 
[21, 22], and a similar result was calculated numerically 
for the p x -p x junction in Ref. [18, 19]. Notice that Eq. 
(16) gives the Josephson current I p (<f>) that is a 2tt- 
periodic functions of <j>, both for T = and T^0. This 
is a consequence of the thermodynamic equilibrium as- 
sumption. At T = 0, this assumption implies that the 
subgap state with the lower energy is occupied, and the 
one with the higher energy is empty. As one can see in 
Fig. 1, the lower energy is always a 27r-periodic functions 



The calculations of the previous section apply in the 
static case, where a given phase difference 4> is maintained 
for an infinitely long time, so the occupation numbers of 
the subgap states have enough time to relax to thermo- 
dynamic equilibrium. Now let us consider the opposite, 
dynamical limit. Suppose a small voltage eV <C Ao is 
applied to the junction, so the phase difference acquires 
dependence on time t: <j)(t) = 2eVt/~h. In this case, the 
state of the system is determined dynamically starting 
from the initial conditions. Let us consider the p x -p x 
junction at T = in the initial state <fi = 0, where the 
two subgap states (11) with the energies ±E are, cor- 
respondingly, occupied and empty. If <p(t) changes suf- 
ficiently slowly (adiabatically), the occupation numbers 
of the subgap states do not change. In other words, the 
states shown by the solid and dotted lines in Fig. 1(c) 
remains, correspondingly, occupied and empty. The oc- 
cupied state (11) produces the current (12): 

T VDeA . /0(t)\ VDeA . ( eVt\ 

(17) 

The frequency of the ac current (17) is eV/h, a half of the 
conventional Josephson frequency 2eV/h. The fractional 
frequency can be traced to the fact that the energies Eq. 
(11) and the corresponding wave functions have the pe- 
riod 47r in (j), rather than conventional 27r. Although at 
4> = 2tt the spectrum in the left panel of Fig. 1 (c) is the 
same as at = 0, the occupation numbers are different: 
The lower state is empty and the upper state is occupied. 
Only at <fi = 4ir the occupation numbers are the same as 
at 4> = 0. 

The 4-7T periodicity is the consequence of the energy 
levels crossing at <f> — it. (In contrast, in the s-wave 
case, the levels repel at <j> = it in Fig. 1(b), thus the en- 
ergy curves are 27r-periodic.) As discussed at the end of 
Sec. IV A, there is no matrix element between the cross- 
ing energy levels at <fi = ir. Thus, there are no transi- 
tions between them, so the occupation numbers of the 
solid and dotted curves in Fig. 1(c) are preserved. In 
order to show this more formally, we can write a gen- 
eral solution of the time-dependent BdG equation as a 
superposition of the two subgap states with the time- 
dependent <f>(t): ip(t) = J2a C a( t )' l PaWt)}- The matrix 
element of transitions between the states is proportional 
to (^(V'+I^V-) = Mi>+\d^H\^)/(E_ -E+). We found 
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that it is zero in the p^-wave case, thus there are no tran- 
sitions, and the initial occupation numbers of the subgap 
states at (f> — are preserved dynamically. 

As one can see in Fig. 1(c), the system is not in the 
ground state when n < <f> < 37r, because the upper energy 
level is occupied and the lower one is empty. In princi- 
ple, the system might be able to relax to the ground state 
by emitting a phonon or a photon. At present time, we 
do not have an explicit estimate for such inelastic re- 
laxation time, but we expect that it is quite long. (The 
other papers [18, 19, 21, 22] that assume thermodynamic 
equilibrium for each value of the phase </> do not have 
an estimate of the relaxation time either.) To observe 
the predicted ac Joscphson effect with the fractional fre- 
quency eV/h, the period of Josephson oscillations should 
be set shorter than the inelastic relaxation time, but not 
too short, so that the time evolution of the BdG equation 
can be treated adiabatically. Controlled nonequilibrium 
population of the upper Andreev bound state was re- 
cently achieved experimentally in an s-wave Josephson 
junction in Ref. [32]. 

Eq. (17) can be generalized to the case where initially 
the two subgap states are populated thermally at (f> = 0, 
and these occupation numbers are preserved by dynami- 
cal evolution 

Ip{f) = | E «H /[W = 0)] (18) 

= sinf^^tanhf^V (19) 
V h J ^DeR \ 2T J V ' 

Notice that the periodicities of the dynamical equation 
(19) and the thermodynamic Eq. (16) arc different. The 
latter equation assumes that the occupation numbers of 
the subgap states are in instantaneous thermal equilib- 
rium for each <\>. 

D. Tunneling Hamiltonian approach 

In the infinite barrier limit D — > 0, the energies ±£q , 
of the two subgap states (11) degenerate to zero, i.e. they 
become midgap states. The wave functions (8) simplify 
as follows: 

<A ±0 = ^° (x)T / fl ° (x) , (20) 
v2 

V>lo = V2k sin(fc F a;) e KX Q ^ 6{-x), (21) 

i> R0 = V2^sin(k F x)e- KX ( _f,_ i0/2 J 9{x). (22) 

Since at D = the Josephson junction consists of two 
semi-infinite uncoupled p^-wave superconductors, ipLo 
and ipuo are the wave functions of the surface midgap 
states [15] belonging to the left and right superconduc- 
tors. Let us examine the properties of the midgap states 
in more detail. 



If (u, v) is an eigenvector of Eq. (3) with an eigenvalue 
E n , then (—v*, u*) for s-wave and (v*,u*) for p-wave are 
the eigenvectors with the energy E n = —E n . It follows 
from these relations and Eq. (1) that j nsky — Ctf n(jk 
with |C| = 1. Notice that in the s-wave case, because 
(u,v) and (— v*,u*) are orthogonal for any u and v, the 
states n and n are always different. However, in the p- 
wave case, the vectors (w, v) and (v* , u*) may be propor- 
tional, in which case they describe the same state with 
E = 0. The states (21) and (22) indeed have this prop- 
erty: 

vlo = iu* L0 , v R0 = -iu* m . (23) 

Substituting Eq. (23) into Eq. (1), we find the Bogoliubov 
operators of the left and right midgap states 

^LOaky = ilLOSky ' ^ROoky = -ilFtnsky ■ ( 2 4) 

Operators (24) correspond to the Majorana fermions dis- 
cussed in Ref. [24]. In the presence of a midgap state, 
the sum over n in Eq. (2) should be understood as 
Sn>o "KV^) Sn=0' wnere we identify the second term 
as the projection Vc of the electron operator onto the 
midgap state. Using Eqs. (23), (24), and (2), we find 

Vc aky {x) = u (x)% aky = v*{x)^l_- k ^. (25) 

Let us consider two semi-infinite p^-wave superconduc- 
tors on a ID lattice with the spacing I, one occupying 
x < I — —I and another x > I. They are coupled by the 
tunneling matrix element t between the sites I and I: 

Hr =T^[cL fey (0c flCTfey (0 + c Raky (l)c Laky (l)]. (26) 

aky 

In the absence of coupling (r = 0), the subgap wave func- 
tions of each superconductor are given by Eqs. (21) and 
(22). Using Eqs. (25), (23), (21), and (22), the tunneling 
Hamiltonian projected onto the basis of midgap states is 

VH T = t [u* L0 (l)u R0 (l) + c.c] (7[ OT 7floT + H.c.) 
= AoVD cos(0/2) (7lo T 7iTO T + 7«) T 7lot), (27) 

where \J~D = 4r sin 2 kpl/hvp is the transmission ampli- 
tude, and we omitted summation over the diagonal index 
k y . Notice that Eq. (27) is 47r-periodic in <j> [24]. 

Hamiltonian (27) operates between the two degener- 
ate states of the system related by annihilation of the 
Bogoliubov quasiparticle in the right midgap state and 
its creation in the left midgap state. In this basis, Hamil- 
tonian (27) can be written as a 2 x 2 matrix 

VH T = A Vd cos(0/2) ( J J ) ■ (28) 

The eigenvectors of Hamiltonian (28) are (1, =f1), i.e. the 
antisymmetric and symmetric combinations of the right 
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and left midgap states given in Eq. (20). Their eigenen- 
ergies are E±{4>) = ^fAoVD cos(^/2), in agreement with 
Eq. (11). The tunneling current operator is obtained by 
differentiating Eqs. (27) or (28) with respect to <f>. Be- 
cause 4> appears only in the prefactor, the operator struc- 
tures of the current operator and the Hamiltonian are the 
same, so they are diagonal in the same basis. Thus, the 
energy eigenstates are simultaneously the eigenstates of 
the current operator with the eigenvalues 

I ± =±— r - sm^-j, (29) 

in agreement with Eq. (17). The same basis (l,Tl) di- 
agonalizes Hamiltonian (28) even when a voltage V is 
applied and the phase <fr is time-dependent. Then the 
initially populated eigenstate with the lower energy pro- 
duces the current I p = \/D(eAo/h) sm(eVt/K) with the 
fractional Josephson frequency eV/h, in agreement with 
Eq. (17). 



E. Josephson current carried by single electrons, 
rather than Cooper pairs 

In the tunneling limit, the transmission coefficient D 
is proportional to the square of the electron tunneling 
amplitude r: D oc r 2 . Eqs. (17) and (29) show that the 
Josephson current in the p x -p x junction is proportional 
to the first power of the electron tunneling amplitude r. 
This is in contrast to the s-s junction, where the Joseph- 
son current (14) is proportional to r 2 . This difference 
results in the big ratio I p /I a = 2/VD between the crit- 
ical currents at T = in the p x - and s-wave cases, as 
shown in Fig. 2 and discussed in Sec. IV B. The reason 
for the different powers of r is the following. In the p x - 
wave case, the transfer of just one electron between the 
degenerate left and right midgap states is a real (non- 
virtual) process. Thus, the eigenenergies are determined 
from the secular equation (28) already in the first order 
of t. In the s-wave case, there are no midgap states, so 
the transferred electron is taken from below the gap and 
placed above the gap, at the energy cost 2Aq. Thus, the 
transfer of a single electron is a virtual (not real) process. 
It must be followed by the transfer of another electron, so 
that the pair of electrons is absorbed into the condensate. 
This gives the current proportional to r 2 . 

This picture implies that the Josephson supercurrent 
across the interface is carried by single electrons in the p x - 
p x junction and by Cooper pairs in the s-s junction. Be- 
cause the single-electron charge e is a half of the Cooper- 
pair charge 2e, the frequency of the ac Josephson effect 
in the p x -p x junction is eV/h, a half of the conventional 
Josephson frequency 2eV/h for the s-s junction. These 
conclusions also apply to a junction between two cuprate 
d-w&ve superconductors in such orientation that both 
sides of the junction have surface midgap states, e.g. to 
the 45°/45° junction (see Sec. V). 




FIG. 3: Schematic drawing of the 45°/45° junction (panel 
a) and 0°/0° junction (panel b) between two d-wave super- 
conductors. The thick line represents the junction line. The 
circles illustrate the Fermi surfaces, where positive and neg- 
ative pairing potentials A are shown by the solid and dotted 
lines. The points A, B, C, and D in the momentum space are 
connected by transmission and reflection from the barrier. 



In both the p x -p x and s-s junctions, electrons trans- 
ferred across the interface are taken away into the bulk by 
the supercurrent of Cooper pairs. In the case of the p x -p x 
junction, a single transferred electron occupies a midgap 
state until another electron gets transferred. Then the 
pair of electrons becomes absorbed into the bulk conden- 
sate, the midgap state returns to the original configura- 
tion, and the cycle repeats. In the case of the s-s junc- 
tion, two electrons are simultaneously transferred across 
the interface and become absorbed into the condensate. 
Clearly, electric charge is transferred across the interface 
by single electrons at the rate proportional to r in the 
first case and by Cooper pairs at the rate proportional 
to t 2 in the second case, but the bulk supercurrent is 
carried by the Cooper pairs in both cases. 



V. JOSEPHSON JUNCTIONS BETWEEN 
d-WAVE SUPERCONDUCTORS 

In this section, we study Josephson junctions between 
two d-wave cuprate superconductors. As before, we se- 
lect the coordinate x perpendicular to the junction line 
and assume that the electron momentum component k y 
parallel to the junction line is a conserved good quan- 
tum number. Then, the 2D problem separates into a 
set of ID solutions (8) in the x direction labeled by 
the index k y . Using an isotropic electron energy disper- 
sion law e = ^(k 2 x + ky)/2m — fj,, we replace the Fermi 
momentum &p and velocity vp by their x-components 

= ^Jk F — ky and vf x — hkF X /m. Thus, the trans- 
mission coefficient D in Eq. (7) becomes fc y -dependent. 
The total Josephson current is given by a sum over all 
occupied subgap states labeled by k y . 

For the cuprates, let us consider a junction parallel to 
the [1, 1] crystal direction in the a-b plane and select the 
x axis along the diagonal [1,1], as shown in Fig. 3a. In 
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these coordinates, the d-wave pairing potential is 

A aky (x, k) = cr2A,3 k v k x /k 2 F , (30) 

where the same notation as in Eq. (4) is used. Direct 
comparison of Eqs. (30) and (4) demonstrates that the 
d-wave superconductor with the 45° junction maps to 
the p x -wave superconductor by the substitution A — > 
cr2Aoky/kF- Thus, the results obtained in Sec. IV for 
the p x -Px junction apply to the 45°/45° junction between 
two d-wave superconductors with the appropriate inte- 
gration over k y . The energies and the wave functions of 
the subgap Andreev states in the 45°/45° junction are 
47r-periodic, as in Eqs. (11). Thus the ac Josephson cur- 
rent has the fractional frequency eV/h, as in Eq. (17). 

On the other hand, if the junction is parallel to 
the [0, 1] crystal direction, as shown in Fig. 3b, then 
&<jk v (x, k x ) = aAp(k% - k 2 )/k F . This pairing poten- 
tial is an even function of k x , thus it is analogous to the 
s-wave pairing potential in Eq. (4). Thus, the 0°/0° junc- 
tion between two d-wave superconductors is analogous to 
the s-s junction. It should exhibit the conventional 2tt- 
periodic Josephson effect with the frequency 2eV/h. 

For a generic orientation of the junction line, the d- 
wave pairing potential acts like p^-wave for some mo- 
menta k y and like s-wave for other k y . Thus, the total 
Josephson current is a sum of the unconventional and 
conventional terms [20] : 

/ = Ci sin(^)/2) + C 2 sin(0) + . . . , (31) 

with some coefficients C\ and Ci- We expect that both 
terms in Eq. (31) are present for any real junction be- 
tween d-wave superconductors because of imperfections. 
However, the ratio C1/C2 should be maximal for the 
junction shown in Fig. 3a and minimal for the junction 
shown in Fig. 3b. 

VI. EXPERIMENTAL OBSERVATION OF THE 
FRACTIONAL AC JOSEPHSON EFFECT 

Conceptually, the setup for experimental observation 
of the fractional ac Josephson effect is straightforward. 
One should apply a dc voltage V to the junction and 
measure frequency spectrum of microwave radiation from 
the junction, expecting to detect a peak at the fractional 
frequency eV/Ti. To observe the fractional ac Joseph- 
son effect predicted in this paper, it is necessary to per- 
form this experiment with the 45°/45° cuprate junctions 
shown in Fig. 3(a). For control purposes, it is also desir- 
able to measure frequency spectrum for the 0°/0° junc- 
tion shown in Fig. 3(b), where a peak at the frequency 
eV/Ti should be minimal. It should be absent completely 
in a conventional s-s junction, unless the junction enters 
a chaotic regime with period doubling [33]. The high-T c 
junctions of the required geometry can be manufactured 



using the step-edge technique. Bicrystal junctions are 
not appropriate, because the crystal axes a and b of the 
two superconductors are rotated relative to each other in 
such junctions. As shown in Fig. 3(a), we need the junc- 
tion where the crystal axes of the two superconductors 
have the same orientation. Unfortunately, attempts to 
manufacture Josephson junctions from the Q1D organic 
superconductors (TMTSF) 2 X failed thus far. 

The most common way of studying the ac Josephson ef- 
fect is observation of the Shapiro steps [34]. In this setup, 
the Josephson junction is irradiated by microwaves with 
the frequency ui, and steps in dc current are detected 
at the dc voltages V n — nhuj/2e. Unfortunately, this 
method is not very useful to study the effect that we pre- 
dict. Indeed, our results are effectively obtained by the 
substitution 2e — > e. Thus, we expect to see the Shapiro 
steps at the voltages V m = mhto/e — 2mfuv/2e, i.e. we 
expect to see only even Shapiro steps. However, when 
both terms are present in Eq. (31), they produce both 
even and odd Shapiro steps, so it would be difficult to dif- 
ferentiate the novel effect from the conventional Shapiro 
effect. Notice also that the so-called fractional Shapiro 
steps observed at the voltage V\ii — tiuj/Ae correspond- 
ing to n = 1/2 have nothing to do with the effect that 
we propose. They originate from the higher harmonics 
in the current-phase relation / oc sin(2</>). The fractional 
Shapiro steps have been observed in cuprates [35], but 
also in conventional s-wave superconductors [36]. An- 
other method of measuring the current-phase relation in 
cuprates was employed in Ref. [37], but connection with 
our theoretical results is not clear at the moment. 



VII. CONCLUSIONS 

In this paper, we study suitably oriented p x -p x or d-d 
Josephson junctions, where the superconductors on both 
sides of the junction originally have the surface Andreev 
midgap states. In such junctions, the Josephson cur- 
rent /, carried by the hybridized subgap Andreev bound 
states, is a 47r-periodic function of the phase difference 
4>: I cx sin(</>/2), in agreement with Ref. [24]. Thus, the 
ac Josephson current should exhibit the fractional fre- 
quency eV/h, a half of the conventional Josephson fre- 
quency 2eV/h. In the tunneling limit, the Josephson 
current is proportional to the first power of the electron 
tunneling amplitude, not the square as in the conven- 
tional case [21-23]. Thus, the Josephson current in the 
considered case is carried by single electrons with charge 
e, rather than by Copper pairs with charge 2e. The frac- 
tional ac Josephson effect can be observed experimentally 
by measuring frequency spectrum of microwave radiation 
from the junction and detecting a peak at eV/fi. 
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